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Abstract 

Using an approximation scheme to deal with the centrifugal (pseudo-centrifugal) term, we solve 
the Dirac equation with the screened Coulomb (Yukawa) potential for any arbitrary spin-orbit 
quantum number k. Based on the spin and pseudospin symmetry, analytic bound state energy 
spectrum formulas and their corresponding upper- and lower-spinor components of two Dirac par- 
tides are obtained using a shortcut of the Nikiforov-Uvarov method. We find a wide range of 

permissible values for the spin symmetry constant C s from the valence energy spectrum of particle 

^ . 

and also for pseudospin symmetry constant C ps from the hole energy spectrum of antiparticle. 

Further, we show that the present potential interaction becomes less (more) attractive for a long 

(short) range screening parameter a. To remove the degeneracies in energy levels we consider the 
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spin and pseudospin solution of Dirac equation for Yukawa potential plus a centrifugal-like term. 
A few special cases such as the exact spin (pseudospin) symmetry Dirac- Yukawa, the Yukawa plus 



CO 

centrifugal-like potentials, the limit when a becomes zero (Coulomb potential field) and the non- 

relativistic limit of our solution are studied. The nonrelativistic solutions are compared with those 

obtained by other methods. 
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I. INTRODUCTION 

The screened Coulomb (Yukawa) potential is widely used in physics, being a good ap- 
proximation to short-range interactions between charged particles in various areas of physics 
[1,2]. In plasma physics it is known as the Debye-Hiickel potential describes the shielding 
effect of ions embedded in plasmas [3]. It has also been used to play a fundamental role in 
(dusty/complex) plasma and colloidal suspensions. The momentum transfer in pair colli- 
sions of particles interacting via the Yukawa potential is well investigated in the limit when 
the interaction is "weak" in the sense that its range (distance at which the interaction energy 
is equal to the kinetic energy) is much shorter than the plasma screening length. This limit 
is known as the theory of Coulomb scattering and is extensively used to describe collisions 
in usual electron-ion plasma [4] . In solid state, atomic and molecular physics it is called the 
Thomas- Fermi or screened Coulomb potential due to the cloud of electronic charges around 
the nucleus [5,6]. Also this potential is well known in nuclear physics as the dominant central 
part of neutrons-protons nuclear interaction due to the massive field exchange (one pion) 
whose mass is m [7,8] . In high energy physics, the potential is used to model the interaction 
of hadrons in short range gauge theories where coupling is mediated by the exchange of a 
massive scalar meson [1,9]. It is defined as follows [10]: 

V(r) = --e- ar , A>0 (1) 

r 

where a and A are the screening (range) and coupling strength parameters, respectively. The 
two parameters are given by different expressions depending on the type of the problem under 
consideration. For example, A = g 2 is positive for attraction, g denotes the coupling constant 
between meson field and the fermion field with which it interacts. Since the field mediator 
is massive, the corresponding force has a certain range, which is inversely proportional to 
the mass, a = mc/h If the mass is zero, then the Yukawa potential becomes equal to 
a Coulomb potential and the range is said to be infinite. Further, the number of bound 
states of the Yukawa potential is found to be finite. Unfortunately, since the Schrodinger 
equation for the screened Coulomb potential does not admit an exact analytical solution [11], 
therefore, various numerical [12,13] and analytical [14-17] methods have been developed in 
the past. Also the energy spectrum can be calculated with high accuracy by means of the 
hypervirial relations and Pade approximation methods [18,19]. The short-range behavior of 
the decaying exponential factor e~ ar and singularity at r = make the task of obtaining 



accurate solutions a difficult task. Besides, most of these calculations suffer from limited 
accuracy when a wider range of potential parameters are being considered [20] . 

An approximate perturbative method has been developed to obtain the energy spec- 
trum and wave functions of the Schrodinger with the Yukawa-like potentials [21]. This 
method has been applied to obtain the energy spectrum and wave functions for the more 
general exponential-cosine-screened Coulomb potentials. These potentials are containing 
an additional Coulomb term superposed with the Yukawa potential that might be useful 
in describing the effective interaction in many-body problem [22]. Further, the asymptotic 
iteration method is used to obtain the energy eigenvalues of the Yukawa potential [23]. The 
J-matrix method has been applied to the Yukawa potential with no special treatment of 
its singularity by using the oscillator basis and the reference Hamiltonian contained only 
the kinetic energy operator [24] . The bound states spectrum and resonance energies of the 
Yukawa potential have been studied using the method of complex scaling [25]. Therefore, 
an alternative Laguerre basis has also been used to comply with the J-matrix requirement 
of a tridiagonal matrix representation of the reference Hamiltonian that includes the r _1 
singularity [18]. 

When a particle is in a strong interaction (range of interaction exceeds the screening 
length, A = a" 1 ), the relativistic effect must be considered which gives the correction for 
nonrelativistic quantum mechanics. The relativistic treatment is of much interest especially 
when (at least) one of the particles is highly charged and their relative velocity is small. 
Taking the relativistic effects into account, a particle in a potential field should be described 
with the Dirac equation. Therefore, the solution of the Dirac equation can be important in 
different fields of physics like nuclear and molecular physics [7,26] . Within the framework 
of the Dirac equation the spin symmetry arises if the magnitude of the spherical attractive 
scalar potential S(r) and repulsive vector potential are nearly equal (i.e., S(r) ~ V(r)) in the 
nuclei (i.e., when the difference potential A(r) = V(r) — S(r) = C s , with C s is an arbitrary 
constant); however, the pseudospin symmetry occurs if S(r) ~ —V(r) are nearly equal (i.e., 
when the sum potential S(r) = V(r) + S(r) = C ps , with C ps is an arbitrary constant) [27]. 
The spin symmetry is relevant for mesons [28]. The pseudospin symmetry concept has been 
applied to many systems in nuclear physics and related areas [27-31] and used to explain 
features of deformed nuclei [32], the super-deformation [33] and to establish an effective 
nuclear shell-model scheme [29,30,34]. Recently, the spin and pseudospin symmetries have 



been widely applied on several physical potentials by many authors (cf. [35-37] and references 
therein). Many authors have investigated approximately the solution of the Dirac equation 
with a few potential models such as the the generalized Morse potential [35], the Hulthen 
potential [36], the Rosen- Morse potential [37] and the screened Coulomb potential [38] etc 
within the framework of various methods. 

In the framework of the spin symmetry S(r) ~ V(r) and pseudospin symmetry S(r) ~ 
—V(r), the bound state energy eigenvalues and associated upper- and lower-spinor wave 
functions are investigated by means of the Nikiforov-Uvarov (NU) method [39]. We have 
approximately solved the Dirac equation for the Hulthen potential [36] with spin and pseu- 
dospin symmetry for any spin-orbit K state and found the eigenvalue equation and corre- 
sponding two-component spinors within the framework of an approximation to the term 
proportional to 1/r 2 . We have also solved the (3+1) dimensional Dirac equation for a single 
particle trapped in the spherically symmetric generalized WS potential under the conditions 
of exact spin and pseudospin symmetry combined with approximation for the spin-orbit 
centrifugal (pseudo-centrifugal) term, and calculated the two-component spinor wave func- 
tions and the energy eigenvalues for any arbitrary spin-orbit k bound states [40]. Recently, 
Setare and Haidari [41] have solved the Dirac- Yukawa problem in the presence of the spin 
symmetry and given only analytical expressions for energy eigenvalues and wave functions. 
However, they have not given further numerical discussions for the validity of their analytical 
solutions. On the other hand, the subject of the pseudospin symmetry of the Dirac- Yukawa 
problem introduced by Ginocchio [27] has not been investigated by Ref. [41]. Over the 
past years, the interest in the quality of the pseudospin symmetry has been increased in the 
framework of the single-particle relativistic potential models. Therefore, we have found that 
it is necessary to give a detailed study for the solution of the Dirac equation with screened 
Coulomb (Yukawa) potential model in the presence of spin as well as pseudospin symmetry 
in a very simple and elegant way by using a shortcut procedures for the NU method. We 
also give a detailed discussion for the validity of the present numerical as well as analytical 
solutions. We also try to explore the exact relativistic energy spectrum of the Coulombic 
field (when a — 0, the low screening range of the Yukawa potential) under the exact spin 
and pseudospin symmetry. 

The analytic solution of the Dirac equation with the screened Coulomb potential is diffi- 
cult to find due to the centrifugal (pseudo centrifugal) term k(k + l)r~ 2 (k(k — l)r~ 2 ) and 



the singular interactions like r _1 (e.g., the Coulomb potential). Nevertheless, employing the 
approximation provided by Greene and Aldrich [42] to the centrifugal term r~ 2 and to the 
singular Coulombic part r _1 makes the solution handy. We work within the framework of 
the low screening parameter throughout the paper. We find analytically approximate bound 
state solutions including the energy spectra and the corresponding spinor wave functions in 
the presence of the spin symmetry and pseudospin symmetry concept for any /t-state within 
the parametric generalization of the NU method [43] given in Appendix A 

This paper is organized as follows. In section 2, we investigate the bound state energy 
equation and the corresponding two-component spinor wave functions in the presence of 
spin and pseudospin symmetry concept for the screened Coulomb potential by employing 
a parametric generalization of the NU method. In section 3, we study some special cases 
like Schrodinger- Yukawa, Dirac-Coulomb, exact spin (pseudospin) symmetric Dirac- Yukawa 
problem and Yukawa plus centrifugal-like potentials. In section 4, we present some numer- 
ical results to the non-relativistic and relativistic numerical energy levels for the Yukawa 
potential. The relevant conclusions are given in section 5. 

II. THEORETICAL FRAMEWORK OF DIRAC BOUND STATE SOLUTIONS 

The Dirac equation for fermionic massive spin-1/2 particles moving in an attractive scalar 
potential S(r) and a repulsive vector potential V(r) is given by [7] 

[col ■ p + (3 (Mc 2 + S(r)) + V(r) - E] ^(r) = 0, ^ nK (r) = ^(r, 9, 0), (2) 

where E is the relativistic energy of the system, M is the mass of a particle, p = —ih'V is 
the momentum operator, and at and (3 are 4x4 Dirac matrices, i.e., 



a = * , P = , ai= , a 2 = , cr 3 = , (3) 





where I denotes the 2x2 identity matrix and (Ti are the three-vector Pauli spin matrices. 
In a spherical symmetrical nuclei, the total angular momentum operator of the nuclei J and 
spin-orbit matrix operator K = —(3 (<r • L + I) commute with the Dirac Hamiltonian, where 
L is the orbital angular momentum operator. The spinor wave functions can be classified 
according to the radial quantum number n and the spin-orbit quantum number k and can 



be written using the Pauli-Dirac representation in the following forms: 

,„ w = ('~<^=i( F - W5 p> ), ( 4) 

where the upper- and lower-spinor components F nK (r) and G nK (r) are the real square-integral 
radial wave functions, Y^ mK {r) and 7^, Jf] are the spin spherical harmonic functions 
coupled to the total angular momentum j and it's projection m on the z axis and k(k + 1) = 
1(1 + 1) and K (k — 1) = 1(1 + 1). The quantum number k is related to the quantum numbers 
for spin symmetry / and pseudospin symmetry / as 

{-(/ + 1) = - (j + \) , (si /2 ,p 3 /2, etc.), j = l + \, aligned spin (k < 0) , 
(o) 
+/ = + (j + |) , (pi/2, 4/2, etc.), 3 = l-\, unaligned spin (k > 0) , 

and the quasi-degenerate doublet structure can be expressed in terms of a pseudospin angular 
momentum ? = 1/2 and pseudo-orbital angular momentum I which is defined as 

-T= - (j + \) , (si/2, p 3 /2, etc.), j = I- 1/2, aligned spin (« < 0) , 

(T+ l) = + (j + |) , (4/2, /5/2, etc.), j = 7+ 1/2, unaligned spin (« > 0) , 

(6) 
where k = ±1, ±2, • • • . For example, (lsi/2, Oc/3/2) and (2p 3 / 2 , 1/5/2) can be considered as 
pseudospin doublets. 

Upon direct substitution of Eq. (4) into Eq. (2), the two radial coupled Dirac equations 
for the two spinor components can be obtained as 

{tt + ?) FUr) = t M ° 2 + EnK ~ A(r) ] GnK(r) ' (7a) 

j- r -^) Gnn(r) = [Mc 2 - E nK + £(r)] F nK (r), (7b) 

where A(r) = V(r) — S(r) and S(r) = V(r) + S'(r) are the difference and sum potentials, 
respectively. 

In the presence of the spin symmetry ( i.e., A(r) = C s = constant), one can eliminate 
G nK (r) in Eq. (7a), with the aid of Eq. (7b), to obtain a second-order differential equation 
for the upper-spinor component: 



d 2 k(k + 1) 1 /^ ,, 9 _ x _,, 



F„Jr) 
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h 2 c 2 



[E 2 nK - M 2 c A - C s (E nK - Mc 2 )] F n 



and the lower-spinor component is obtained from Eq. (7a): 

1 ( d 



G„.Jr) 



C\ ( lr + ^ F - [n 



(8) 



(9) 



Mc 2 + E. 

where E nK ^ —Mc 2 , only real positive energy spectrum exist when C s = (exact spin 
symmetry). On the other hand, in the presence of the pseudospin symmetry ( i.e., S(r) = 
C ps = constant), one can eliminate F nK (r) in Eq. (7b), with the aid of Eq. (7a), to obtain 
a second-order differential equation for the lower-spinor component: 
72 „(„ r\ i 



d k(k — V, 
+ 



dr 



h 2 c 2 



(Mc 2 - E nK + C ps ) A( 



G n Jr) 



h 2 c 2 



[E 2 nK - M 2 c 4 - C ps (E nK + Mc 2 )] G ri 



and the upper-spinor component F nK (r) can be obtained from Eq. (7b) as 



F n Jr) 



Mc 2 - E, n 



H 



C ps \dr r 



G nK (r), 



(10) 



(11) 



where E nn ^ Mc 2 , only real negative energy spectrum exist when C ps = (exact pseudospin 
symmetry). The physical solution demands that the upper and lower radial components 
should satisfy the boundary conditions: F nK (0) = G nK (0) = and F nK (oo) = G nK (oo) = 0. 

A. Spin symmetry Dirac-Yukawa problem 



At first, we investigate the spin symmetry in the form of £77(2) by taking the S(r) = 
2V(r) — > Vy{r) [44] which can be easily reduced into the non-relativistic limit under a 
certain appropriate transformations. Equation (8) shows that the energy eigenvalues, E nK 
is mainly dependent on the quantum numbers n and /. For example, when I ^ 0, the states 
with j — I ± 1/2 are degenerate. The sum potential S(r) in Eq. (8) is simply taken as the 
Yukawa potential, 



Er 



A>0, 



(12) 



which provides a simple Schrodinger-like equation in the form: 
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F nK (r) = 0, k(k + 1) = 1(1 + 1) (13) 



where 

V \ = A 9I2 9 ( M ° 2 ~ E ™) (MC 2 + E nK - C a ) , Wi = — ^ [#»* + Mc ' ~ C-] A ( 14 ) 

4a z n c l 2an c z 

with k values are given in Eq. (5). The exact analytic solution of Eq. (13) is difficult to 

find due to the centrifugal kinetic energy term k(k + l)r -2 and the singularity of r -1 -type. 

Nonetheless, if k is not too large, the case of the vibrations of small amplitude about the 

minimum, we attempt to use the Greene-Aldrich [42] conventional approximation to deal 

with centrifugal term, 

— ^4a 2 *. (15) 

r 2 (l _ e - 2ar ) 

Introducing the new parameter, x(r) = e~ 2ar G [0, 1] and further substituting Eq. (15) into 

Eq. (13), we obtain 

d 2 (1 — x) d 

dx 2 x(l — x) dx 

+ bM±^)f!±P^^±M£ria \ FUx) = 0, (16) 

X z (l — Xf J 

where we have inserted F nK (r) = F nK (x). In order to clarify the parametric generalization 
of the NU method [35,37,43], let us take the following general form of a Schrodinger-like 
equation written for any potential, 



d 2 t(x) d <j(x) 



il> n (z) = 0, (17) 



dx 2 a(x) dx a 2 (x) 
satisfying the wave functions 

tW = <j)(x)y n (x). (18) 

In addition, the two polynomials 

t(x) = ci - c 2 x, (19) 

and 

<j(x) — x (1 — c^x) and a(x) = — £ 2 :r 2 + ^ ± x — £ , (20) 

are at most of first- and second-degree, respectively. Comparing Eq. (17) with its counter- 
part Eq. (16), we obtain values for the constants q (i = 1,2,3) along with £■ (j = 0, 1,2). 
Now, following the NU method [39] and making the substitution of Eqs. (19) and (20) leads 
to more general forms for the polynomials ff(z) and r(z), the root k, the eigenvalues equa- 
tion and the wave functions (f>(z) and y n (z) all expressed in terms of the constants q (i = 4, 
5, • • • , 13) as given in Appendix A. Hence, the task of computing the energy eigenvalues 



and the corresponding wave functions of Eq. (13) within the framework of the parametric 
generalization of the NU method becomes relatively easy and straightforward. It may be 
explained shortly in the following steps: 

Firstly, we need to find the specific values for the parametric constants q (i — 4, 5, • • • , 13) 
by means of the relation Al of Appendix A. The values of all these constants q (i — 1, 
2, • • • ,13) together with £• (j = 0,1,2) are therefore listed in Table 1 for the screened 
Coulomb potential model. 

Secondly, by using the relations (A2-A5), the analytic forms of the essential polynomials 
n(x) and t(x) along with the root k, required by the NU method, can also be found as 



n(x) = -\~\ [(2« + 1 + 2i/x) x - 1v x \ 



k = —k (k + 1) + u>i — (2k, + 1) Ui, 



and 



t(x) = 1 + 2i/i - 2 ( 1 + v x + - (2k + 1) ) x, 



(21) 
(22) 

(23) 



where r'(x) < must be satisfied in order to obtain physical solution according to the NU 
method [39]. 

Thirdly, we need to calculate the energy eigenvalues by means of the eigenvalue equation, 
relation A6 which gives 



(n + k + 1) + 2 (n + k + 1) v x = ui. 



(24) 



Finally, after making use of Eq. (14), the above equation for the Yukawa potential can be 
expressed implicitly in terms of the energy E nK as 

^(Mc* - E nK ) (M? + E nK - C.) + ahc(n + K.+ l)= (Mc ' + f ^ ' 9A^L , 

Mc 2 > E nK and Mc 2 + E nK > C s , n = 0, 1, 2, • • • , (25) 

where k ^ — (n + 1) and A = Ze 2 . The above energy equation can be also rearranged in a 
quadratic form (in relativistic units H — c — 1) as 
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K*i, 
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) w 



E nK + ( -^ J 



aiVi 



(M + aA) W = 0, 
(26) 



where 

iVi = 2 (n + k + 1) , W = C s - M, S = ^(C s + aA). (27) 

The two energy spectrum formula of the quadratic equation (26) is 



(A 2 W + SN 2 ) ± J (A 2 W + SN?) 2 - (A 2 + N 2 ) (AW + \aN 2 ) + MWN 2 

(28) 
where {A 2 W + SN 2 ) > J (A 2 + N 2 ) (AW + \aN 2 ) 2 + MWN 2 for distinct particle and 
anti-particle real energy bound states E? K = E+ K and E® K = E~ K , respectively. Otherwise, 
in the case when the above inequality does not hold, we will have no bound state solutions 
(scattering states). 

Figure 1 shows the ground state valence energy level of particle and antiparticle as a 
function of different values of the coupling constant A and the screening parameter (range) 
of the potential a = 0.01, 0.02, 0.05 and 0.10 for two degenerate partners (n — 0, k — 1) and 
(n — 3, k — —2) labelled as (Opi/2, 3^3/2) of particle and antiparticle with the parameters 
choices of mass M = 5.0 fm~ l and spin constant C s = 4.9 /to -1 , as it should be expected, 
for a given value of a the bound state becomes sharply (slowly) deeper (more attractive) 
for particle (antiparticle) on increasing the coupling constant A (heavy nucleus). On the 
decreasing the value of the range a (lower screening range), the energy goes to a more neg- 
ative value for antiparticle and to a less positive value for particle (energy becomes more 
attractive). On the other hand, for a fixed value of A, the bound state energy becomes shal- 
lower on increasing a for particle. It is seen that in the limit of a very short-ranged potential 
(a — ¥ 0), the potential approaches the S- function limit that can bind particles and antipar- 
ticle stronger than finite-ranged potentials (1). Furthermore, Fig. 2 shows the variation of 
ground state valence energy level of particle and antiparticle with the spin constant C s for 
several values of the spin-orbit k — 1,3 and 5 with special choices of parameters M = 5.0 
/to -1 , a = 0.1 /to -1 and A = 1. A very careful inspection for both numerical results and 
Fig. 2a shows that there is a small energy difference between the states k — 1,3 and 5 
although the values of spin constant C s increases in the range 11.4 /to -1 < C s < 20 /to -1 , 
i.e., the energy spectrum is not sensitive to the influence of k in the aforementioned range. 
With the increasing k value, we see that Eq k fan out toward the stronger positive energy 
spectrum for the given range —20 /to -1 < C s < 8.8 /to -1 . The scattering (not bound) states 
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could be seen in the range 8.8 fm" 1 < C s < 11.3 /m _1 for k — 1,3 and 5. The physical 
explanation to Figure 2a is being illustrated as follows. It is seen in Table 2 that the energy 
spectrum (positive and negative) of Eq. (25) is entirely dependent on the choice of C s , 
i.e., Eq k = Eq k (C s ). For example, in the single electron interaction with nucleus Ze (units 
H — c — e — 1), we may have two essential requirements: (a) M > E nK and M + E nK > C s , 
in which both impose a restriction on the choice of the range values of the permissible spin 
symmetry constant C s are in the interval —20 fm -1 < C s < 8.8 fm" 1 with a requirement 
that energy spectrum be real along with \E nR \ < M and (b) M < E nK and M + E nK < C s 
which is possible for 11.4 fin" 1 < C s < 20 fm" 1 where E nK > M. In the two cases we 
have taken the spin-orbit quantum number k — 1, 3 and 5. While taking k — 1, 3, 5, 7 and 
9, C s lies in the two intervals: -20 fm" 1 < C s < 8.0 fm" 1 and 12.2 fm" 1 < C s < 20 
fm~ l . Nonetheless, it is necessary to choose the physical solution as shown in case (a) since 
it is consistent with the exact spin symmetry when C s = (i.e., S(r) = V(r)) (\E nK \ < M) 
in which the energy levels become very sensitive to the influence of k as usually expected. 
Table 2 supports our choice of the allowed range —20 fm^ 1 < C s < 8.8 fwT 1 . Moreover, 
we find out that the second case is not sensitive to the influence of k. Thus, for our choice of 
C s = 4.9 fm" 1 which is falling in (a), we present all ground energy spectrum including four 
states: E ,i = 4.78641 fm,- 1 , £ ,2 = 4.94209 fm" 1 , £ 0)3 = 4.98904 fm" 1 and E 0A = 4.99998 
fm" 1 . 

The results presented in Fig. 2b and Table 2 show that the energy difference of antiparticle 
between the states k — 1,3 and 5 shows a slight change although the values of spin constant 
C s increases in the interval 10.6 fm" 1 < C s < 20 fm" 1 (same as in Fig. 2a). This 
range is forbidden since Eq k < —M. As K decreasing, Eq k fan out toward the stronger 
(deeper) negative energy spectrum when —16 fm" 1 < C s < —0.1 fm" 1 (Eq k is very strongly 
negative, forbidden). In the intervals fm" 1 < C s < 4.9 fm" 1 , Eq k < M (negative, 
permissible), 5.0 fm" 1 < C s < 9.6 fm" 1 , Eq k < M (positive, permissible) while 9.7 
fm" 1 < C s < 10.6 jm' 1 (energy is complex). The variation of the ground negative energy 
level with C s contains forbidden values of C s that gives Eq )K ^ — M which results in an 
infinity wave function. The range values of the allowed C s are in the interval fm" 1 < 
C s < 9.6 fm" 1 with E 0)K ^ -M when C s = 0. 

Next, in order to establish the wave functions F nK (r) of Eq. (8), the relations (A7-A10) 



11 



are used. Firstly, we calculate the first part of the wave functions, 

(J)( x ) = x "i(i _ x y+\ Vl > o. (29) 

The weight function takes the form 

p(x)=x 2ui (l-x) 2K+1 . (30) 

which can generate the second part of the wave functions, 

d n 



yjx) ~ x~ 2v H\ - x)-^ +1 ^^— \x n+2l/1 (1 - x) n+2K+l ] « p(2"i-2«+i)(i _ 2a: ), (31) 



where Pi' (1 — 2x) is the orthogonal Jacobi polynomials [45,46]. Finally, the upper spinor 
component F nK (x) for any arbitrary k can be obtained by means of Eq. (18) as 



F n Jr) = Af nK e- 2 ^ ar (l - e -^ r+ i p (2u 1: 2 K+ i) (1 _ 2e 



2ar\ 



p( 2 , 1 , 2K+ i) (1 _ 2e - 2ar) = r (" + 2 ^ + )) aFl (-n,n + 2(i/ 1 + K + l);l + 2i/ 1 ; e -*"-), (32) 

1 yZUi ■+" 1)71! 

where the normalization constants A/" nK are calculated in Appendix B. 
The derivative relation of the hypergeometric function, 



d 

dx 



,ab\ 

2 Fi (a; o; c; x) = ( — ) 2 F 1 {a + 1; 6 + 1; c + 1; x) , 



is usually used to calculate the corresponding lower-component G nK {r) by means of Eq. (9): 



M nK (2ot(K + l)e "" , k 



-2ar 



G nK (r) = — i — —i- 2ai/i + - F n/S (r) 

y) Mc 2 + E nK -C s \ l-e~ 2ar r) y! 

,»/■ 2wa (w + 2t/i + 2k + 2) _ 2anis+ i , _ 2a r\"i+i 

+ ' V " K (Mc2 + P nK -C,)(l + 2z. 1 ) 1 j l j 

x 2^1 (-n + 1; n + 2 (i/i + k + 1) + 1; 2 (1 + ^i) ; e~ 2ar ) . (33) 

The hypergeometric series 2 P X (— n + 1; n + 2 (z/i + k + 1) + 1; 2 (1 + z/i) ; e~ ar ) is termi- 
nated for n = and thus does not diverge for all values of real parameters v\ and k + 1. 
In relativistic units ft = c = 1, P nK 7^ —M, i.e., negative energy states are forbidden, when 
C s = in which the positive energy solution of Eq. (25) is required. Therefore, the negative 
solution is not desirable, see Fig. 2b.. 

In Fig. 3, we plot the upper and lower spinor wave functions of ground Opi/2 and first 
excited lpi/ 2 states for (a) particle and (b) antiparticle with k — 1. The choices of parameters 
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M = 5.0 fm -1 , C s = 4.9 fm~ l , a = 0.1 fmT 1 and A = 1 are used. In case of the particle 
(positive energy), the upper and lower spinor wave functions of the ground (first excited) 
state are found to be similar in shape. It is noted that the amplitude of the upper wave 
function F^r) (F^(r)) is nearly three times larger than the lower wave function (joi(r) 
(Gn (?"))■ Further, the amplitude of the ground state wave function F^[(r) is nearly two 
times larger than the first excited state wave function F^[(r). The range of the upper 
component is wider than the lower component of the wave function and the range of the 
first excited state is wider than the range of the ground state wave function. 

Let us now study the nonrelativistic case. Making the appropriate changes: v\ — > e 2 nl = 
- (2mE nl /4a 2 h 2 ) ,u 1 ->'y 1 = (mA/ah 2 ) and k(k + 1) ->■ 1(1 + 1) in Eqs. (16), (21)-(24) and 
(29)-(32) together with Table 1, we can easily obtain the energy spectrum of the Schrodinger 
equation for the Yukawa potential model: 

7i (n + l + 1) 
£nl ~ 2(n + / + l) 2 ' (34) 



which can be explicitly expressed as 



h 



Enl ~ "2^ 



and the radial wave functions: 



mA 



h 2 (n + l + 1) 



(n + l + 1) a 



2 



(35) 



R nl (r) = r- l u nl (r) = N^r' 1 e~ €nir (1 - e -^y+i p (2s nl ,2i+i)^ _ 2e - 2ar ) 
K ' %tr+\|n! ) e ' £ "' r(1 " 6 ' 2ar)l+1 2Fl (_n ' U + 2 {£nl + l + i y> l + 2£ ^ ^^ 



where e nt = 2ae ni 



^TT)-(- + ^ + 1 )« 



(36) 
and the normalization constants M n i are 



carried out in Appendix B. 



B. Pseudospin symmetry Dirac- Yukawa problem 

From Eq. (10), we can see that the energy eigenvalues depend mainly only on n and /, 
i.e., E nK = E(n,l(l + 1)). For I ^ 0, the states with j — I ± 1/2 are degenerate. This is a 
SU(2) pseudospin symmetry. We impose A(r) as the Yukawa potential model given in (1): 

A(r) = -^V Qr , A = A > (37) 

13 



leading to a Schrodinger-like equation in the form: 
dr 2 



r 






G n Jr) 



0, 



(38) 



with 



V l = A 9I2 o ( Mc2 + E ™) ( Mc2 ~ E nn + C ps ) , 002 = ^4^T l E ™ - M ° 2 ~ C Ps] A , (39) 

4a 2 h c 2 2an c 2 

where k(k — 1) = 1(1 + 1) satisfying Eq. (6). We follow the same procedures of solutions 
discussed before to obtain a Dirac equation satisfying G nK (r), 

d 2 (1 — x) d 



dx 2 x(l — x) dx 

+ [-(^ + ^)^ + (2^ + a, 2 -«(«-l))x-^]| = Q 

x 2 (l — x) 2 J 

To avoid repetition in the solution of Eq. (40), a careful inspection for the relationship 

between the present set of parameters (u>2, v\) and the previous one (u±, u\) tells us that the 

negative energy solution for pseudospin symmetry, where S(r) ~ —V(r), can be obtained 

directly from the spin symmetric solution by using the following parameter mapping [40,47]: 

F nK {r) <-► G nK (r), k -> «-l, V(r) ->■ -V(r) (i.e., A ->■ -A), £ nK ->■ -£ nK and C s -»■ -C ps . 

(41) 
Further, the constants in the case of pseudospin symmetry concept are listed in Table 1. 
Applying the above transformations to Eqs. (21)-(24) leading to the following pseudospin 
symmetric energy equation, 

(n + k) 2 + 2(n + K)u 2 = w 2 - (42) 

Finally, with the aid of Eq. (39), Eq. (42) can be also expressed in terms of the energy, 

(E nR -Mc 2 -C ps )A 



(Mc 2 + E nR ) (Mc 2 - E nR + C ps ) + ahc (n + K 



2hc (n + k) 



Mc 2 > -E nK and Mc 2 + C ps > E nK , n = 0, 1, 2, 



(43) 



where k 7^ — n. When C ps = (exact pseudospin symmetry, S(r) = —V(r) case) then we 
require — M < E nK < M.The above pseudospin energy equation can be rearranged in a 
quadratic form (H — c — 1): 



A 



E — 2 



T 






[/ 



EL, 



A£/\ 2 faN 2 x2 



N? 



■(M-aA)U = 0, (44) 
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where 

N 2 = 2 {n + k) , U = C ps + M, T = ^(C ps + aA). (45) 

The above quadratic energy equation can be easily obtained by means of Eq. (25) through 
making the replacements: n + K + l — > n + k,, A — > —A, E nK — > —E nK and C s —¥ —C ps . The 
two energy solutions of the quadratic equation (44) can be obtained as 



(A 2 U + TN 2 ) ± J (A 2 U + TN 2 ) 2 - (A 2 + N 2 ) (AU + \aN 2 ) - MUN, 

rp± _ V [ ' 

A 2 + N 2 



(46) 



where (A 2 U + TN 2 ) > J (A 2 + N 2 ) (AU + \aN 2 ) - MUN 2 for distinct particle ££„ and 
anti-particle E® K real bound state energies. 

On the basis of pseudospin symmetry, Fig. 4 shows the ground state hole energy level 
of particle and antiparticle for different values of the coupling constant A and the screening 
parameter (range) of the potential a = 0.01, 0.02, 0.05 and 0.10 for two degenerate partners 
(n = 0,k = 2) and (n = 1,k — —3) labelled as (0^3/2,1^5/2) with M = 5.0 fm^ 1 and 
C ps = —5.0 jvrC x . In this exploratory investigation, as it should be expected, for a given 
value of a the bound state becomes slowly (sharply) deeper for particle (antiparticle), i.e., 
it becomes more (less) attractive, on increasing the coupling constant A (heavy nucleus). 
On increasing the value of the screening parameter a, the particle (antiparticle) becomes 
more (less) attractive as A increasing. Therefore, in the pseudospin symmetry case, the 
antiparticle is less attractive to heavier nuclei while in the spin symmetry case the particle 
is more attractive to heavier nuclei. In addition, increasing the value of screening parameter 
a in both particle and antiparticle leads to less attractive interaction with heavier nucleus. 

Furthermore, Fig. 5 shows the variation of the ground state hole energy for particle and 
antiparticle as a function of different values of pseudospin constant C ps plotted for several 
values of the spin-orbit k — 1, 3 and 5. We take the set of parameters M = 5.0 /m _1 , a — 0.1 
fm^ 1 and A — 1. The results presented in Fig. 5a show that the energy difference of particle 
between the states k — 1, 3 and 5 is almost same although the values of spin constant C ps 
increases in the interval —8 fm" 1 < C ps < 20 fwT 1 . With an increasing of k value, Eq k 
is seen to fan out toward the stronger negative energy spectrum (more attractive) when 
—20 ]vrT x < C ps < — 13 ]vrT x . The scattering states could be seen in the short interval 
— 12 fm^ 1 < C ps < —8 fwT 1 . Figure 5b shows that the energy difference of antiparticle 
between the states n = 1,3 and 5 is almost same although the values of spin constant C, 



p.s 
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increases in the interval —20 fm~ x < C ps < — 11 fm^ 1 (i.e., antiparticle energy state is 
not sensitive to spin-orbit quantum number k). With the k increasing, Eq k fan out toward 
the stronger (deeper) negative energy spectrum when —8 fm^ 1 < C ps < 20 ]vrT x . The 
scattering states could be seen in the short interval —11 ]vrT x < C ps < —8 fwT 1 . In the 
presence of pseudospin, we consider the physical case where the energy has to be negative 
(i.e., E nK 7^ M when C ps = 0) which is simply the case of antiparticle. As seen in Fig. 5b, 
the range of the allowed pseudospin constant C ps falls in the range —8 fm~ l < C ps < 20 
fm~ l in which E~ K is sensitive to the influence of k and also energy results satisfy the 
condition \E nK \ < M. Consequently, we choose C ps = —5.0 fm^ 1 in the present numerical 
calculations. 

We present the essential procedures in calculating the wave functions. The first part of 
the wave functions is 

(f)(x) =x U2 (l-x) K , u 1 >0 (47) 

and the weight function is 

p(x)=x 2u *(l-x) 2K - 1 . (48) 

and this generates the second part of the wave functions, 

yJx) ~ x- 2v H\ - x)-^- 1 ^ \x n+2 ^ (1 - x) n+2K - 1 ] » p(2"»-2«-i)(i _ 2x ). (49) 

dx n 

Finally, the lower spinor component G nK (x) for arbitrary k can be obtained by means of Eq. 

(18) as 

G nK (r) = M nK e- 2 ^ m \l - e -2°r)"pe„ 2) 2,e-l) (l _ 2e -2ar )5 

p^^-i )(1 _ a.-^ = r^+ 2, 2 +j) 2 ^ ( _ n>n + 2(l/2 + /6);1 + 2l/2;e - 2 «- ) , (50) 

The upper-component E nK (r) can be calculated from Eq. (11) as follows 

F - W = (Mc 2 -E nK + Cps ) (l^^ ~ *** - r) G - (r) 

, »/■ 2na (n + 2z/ 2 + 2/t) _ 2 , _2ar\"2+i 

+ ^ nK (Mc 2 -E nK + C ps )(l + 2u 2 ) [i } [e } 

x 2^1 (-n + l;n + 2(z/ 2 + /t) + l;2(l + z/ 2 );e- 2ar ). (51) 

In relativistic units h = c = 1, E nK ^ M, i.e., positive energy states are forbidden, when 
C ps = in which the negative (antiparticle) energy solution of Eq. (43) is required. 
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Therefore, the positive solution is not desirable, see Fig. 5a. The hypergeometric series 
2-Fi (—n + 1; n + 2 (z/ 2 + k) + 1; 2 (1 + i/ 2 ) ! e~ ar ) terminates for n = and thus does not 
diverge for all values of real parameters v 2 and k. 

In Fig. 6, we plot the upper and lower spinor wave functions of ground 0^/2 and first 
excited 1^3/2 states for (a) particle and (b) antiparticle with k — 2. The set of parameters 
M = 5.0 /to -1 , C ps = —5.0 /to -1 , a = 0.1 /to -1 and A = 1 are used. 

III. A FEW SPECIAL CASES 



Let us study four special cases. At first, we study the nonrelativistic (Schrodinger- 
Yukawa) case by setting C s = (exact spin symmetry limit) and making the changes k —> I, 
M + E nK — > 2fi/h 2 and M — i?„ K — > —E nK . Hence, from Eq. (25), it follows that 

,2 -.2 



E n l — — TW 



A 



h z 



2// 

¥ [2 (n + I + 1) 2/x 



(n + / + 1 ) a 



, n = 0,l,2,..., / = 0,1,2, 



(52) 



Further, in the limit when a — 0, the above result reduces to the well-known spectrum 



(C) 



Mo 






with a wave functions R n i(r 



— — r— : T > 0, and L'Lix) is the Laguarre 

ft (n+Z+l) ra\ / o 



for the nonrelativistic Coulombic field, _E ni 
r-^nzW = r- 1 e^ r r l+1 L ( a l+1) (2l3r) with /3 
function [22]. 

Second, spin symmetry Dirac-Coulomb (a — > 0) case 

^(Mc 2 - E nK ) {Mc 2 + E nK - Cs) = {M f. ± % ~ ^ A , n = 0,l,2,. 

Arte \Ti ~\- K ~r J- ) 

The above energy equation can be rearranged in a quadratic form {% — c — 1): 



(53) 



1 



A 

N~i 



E' 



g. + 2[A] »• 



^BK *T 



AWY 



iV- 



r 



+ Mlf = 0, 



(54) 



and the two energy solutions of the above equation can be obtained as 

n± (2A 2 W + C S N 2 ) ± y/(2A 2 W + C S N 2 ) 2 - 4W{A 2 + N 2 ) [A 2 W + MN 2 



El 



(55) 



2 {A 2 + N 2 ) 

In the limitation of pseudospin symmetry, the Dirac-Ykawa problem reduces to Dirac- 
Coulomb problem when a — >■ 0, 

(£ nK - Mc 2 - C„,) A 



(Mc 2 + E nK ) (Mc 2 - E nK + C } 



V s ) 



psj 



2hc (n + k) 



n = 0,1,2, 



(56) 
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and it can be rearranged in a quadratic form (h = c = 1) as 



1 + l ^ 



E„ 



c;,/r + 2(A] v 



■l-'tik *T 



At/ 



MU = 0. 



Thus, the two energy solutions of the above equation can be readily obtained as 



, _ (2A 2 U + C ps N 2 ) ± y/(2A 2 U + C pfl iVf ) 2 - 4U{A 2 + JV|) [A 2 £7 - MN, 
nK ~ 2 (A 2 + iV|) 



(57) 



(58) 



Third, on the base of the exact spin symmetry (C s = 0), the energy equation for Dirac 
Yukawa problem becomes (in units h — c— 1) 

Ay/M + E n 



M - E nK + a (n + k + 1 



\E nK \ <M, n = 0,1,2, 



2(ra + /c + l; 
with two energy solutions of the quadratic equation (25) can be obtained as 



(59) 



(aN 2 - 2MA)A ± J (aN 2 - 2MA) 2 A 2 - (A 2 + N 2 ) (aN 2 - 2MA) 2 - (2MN ± ) 



El 



2 (A 2 + N 2 ' 



(60) 



where Ni is defined in Eq. (27). From Eqs. (32) and (33), the upper and lower wave 
functions are 



"nny.') J*nn. 



and 



r(n + 2 7 + l) 
T(2 7 + l)n! 



-2yar i 



1 -e 



— 2ar\K+l 



.F, (-n,n + 2( 7 + / t + l);l + 2 7 ;e- 2OT ), 

(61) 



G nK (r) 

+Afn 



Muk (2a(K+l)e 



-2ar 



M + E nK \ 1 - e 
2na (n + 27 + 2/t + 2) 
'■ (M + E nK )(l + 2 1 ) 



-2ar 



- 2a 7 + - F„ K (r 



2ar\K+l / -2or\7+l 



1 - e' 2m ') K+1 (e~ 2ar ) 



x 



2 Fx (-n + 1; n + 2 (7 + « + 1) + 1; 2 (1 + 7) ; e" ar ) , M ^ -E n 



(62) 



respectively, where 7 = ^a/M 2 — i? 2 K . 

In view of the exact pseudospin symmetry (C ps = 0), the energy equation for Dirac- 
Yukawa problem reads (H — c — 1) 



v/M 2 -E 2 K + a(n + K)= ( ^; M) / , n = 0,l,2, 

2(n + k) 



(63) 



where l-E^I < M and with two energy solutions: 



E 



± 



aN 2 + 2MA)A ± \\(aN 2 + 2MA) 2 A 2 - {A 2 + N 2 ) (aN 2 + 2MA) Z - (2MN 

2 (A 2 + N 2 ) 



(64) 
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The lower and upper-component wave functions are 
G nK (r) = JC ^ h±^ e -^r {1 _ e -2ar r ^ (_ n> n + 2 ( 7 + «) ; l + 2 7 ; e" 2 ") , (65) 



and 



Mn 



-2ar 

! \! . . . . / •■' ft f-» f 

FnK ^ y (Mc 2 -E nK ) \l-e~ 2ar 



r2nHt " -2cry-^)6Ur) 



-a/;. 



2na (n + 2 7 + 2k) _ g _ 2 /_ 2c ^7+i 



(Mc 2 -£ nK )(l + 2 7 ) 



2 F X (-n + 1; ra + 2 (7 + «) + 1; 2 (1 + 7) ; e" ar ) . 



(66) 



respectively. 

Fourth, let us find the analytic solution of the Yukawa plus an added centrifugal-like 
term, i.e., V(r) = — e~ ar + %. The aim behind this choice is to remove the degeneracy of 
energy for various states. In view of spin symmetry, after inserting S(r) = — ^-e~ ar + ^ into 
Eq. (8), we obtain 



-t-^ ^-5 --4a 2 u 2 1 + 2au 1 -e ar 

dr z r z r 



F nK (r) = 0, 



(67) 



and, hence, the energy spectrum can be obtained from Eq. (28) by making the change 



Kj 7 Kl dbS 



y/(M - E nK ) (M + E nK -C s )+a(n + n' + l)= ( M + E ™ . J^A , K '^-(n + l), (68) 



where 



1 1 

-2 + 2 



(2k + l) 2 + AD (M + E nK , - C s 



(69) 



Further, the upper and lower wave functions can be obtained simply via Eqs. (32) and (33) 



as 



F nK (r)=Af n 



-2v\a.r 1 



1-e 



-2ar\K,'+l 



s Fi (-n,n + 2(z/ 1 + /t' + l);l + 2z/ 1 ;e- 2ar ), (70) 



and 



G„Jr) 



M n 



M + E nK -C s \ 1-e 



2a(n' + l)e' 2ar 



-j—2ar 



2avi H F n Jr 



, 2na(n + 2u l + 2K , + 2) 

+MrtK (M + E nK -C s )(l + 2v 1 ) [l 



_ -2ar\n'+l I -2ar\ Vl+1 



2 Fi {-n + 1; n + 2 (y x + «' + 1) + 1; 2 (1 + i/x) ; e~ 2ar ) . 



(71) 
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respectively. In view of pseudospin symmetry, we obtain the energy spectrum from Eq. (46) 
by making the change K — > k" as 

(M + E nK ) (M - E nK + C ps ) + a(n + «") = ( E ^ ~ M~ C ) A ^ ^ = _^ ^ 

where 



*" = \ + \ \A 2k " I) 2 + 4£ (E nK -M- C ps ). (73) 

The lower and upper wave functions can be obtained as follows 

G nK (r) = K tK e~ 2 ^ ar (l - e- 2ar ) K " 2 F, (-n, n + 2 {v 2 + «") ; 1 + 2z/ 2 ; e" 2ar ) , 

and 

F "« (r) = (M-^ + c p ,) (l^r " 2-, " 7 J CM 

,»/■ 2na (n + 2t/ 2 + 2k") „-2a n «" / -2^2+1 

+ ^(M- J E nK + C ps )(l + 2z, 2 ) li j ^ j 

x 2 Fi (-n+l;n + 2(z/ 2 + /t") + l;2(l + z/ 2 );e- 2ar ). (74) 

Fifth, the exact Dirac-Coulomb problem (a = 0) has the following energy equations: [1] 

4iV 2 s (M + £ nK ) = (M - K«) A 2 , iV ps = n + k, n = 0, 1, 2, • • • , (75) 

and 

4iV 2 (M - K«) = ^(M + £„*), N s = n + « + 1, ra = 0, 1, 2, • • • , (76) 

in the limitation of the exact pseudospin (C ps = 0) and spin (C s = 0) symmetry, respectively. 
Obviously, in making the following changes E nK — y —E nK and k — y K — 1, one can easily 
switch off from spin symmetry, Eq. (76), into pseudospin symmetry, Eq. (75). Furthermore, 
Eqs. (75) and (76) are identical to Eqs. (37) and (47) of Ref. [48] (if one sets A = = C 
and B — y A). They are also identical to Eqs. (40) and (52) of Ref. [49]. 

IV. NUMERICAL RESULTS 

From Eq. (35), for small potential strength parameter A = v2, we calculate some 
numerical values of the bound state energies of the Schrodinger equation with the Yukawa 



[1] The bound state solutions are exact for the case a = since it lies in the short screening range. 
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potential for various values of quantum numbers n, I and screening parameter a (h — \x — 1). 
Our approximated results in Table 3 are compared with those of [16,17,19] together with 
the results of [13] who solved the Schrodinger equation numerically. Our results are in good 
agreement for small a values and becomes worse as a increases since the approximation 
used to substitute the singular part of Yukawa potential, r _1 and the centrifugal term r~ 2 
are true for or <§; 1. As the potential strength parameter A becomes larger, the numerical 
solution of the Schrodinger equation is extremely difficult because the screened Coulomb 
(Yukawa) potential is very deep and the wavefunction becomes very sharply peaked near 
the origin. Because of the instability of the wave function the energy eigenvalues become 
fairly inaccurate as the strength A increases. 

Based on the spin symmetry case, from Eq. (28), we can calculate some relativistic 
particle E+ K and antiparticle E~ K bound state energies with values of parameters a — 0.1 
/to -1 , A = 1.0, M = 5.0 fin" 1 and C s = 4.9 /m _1 for various states with quantum numbers 
n and I (k < 0, k > 0) in the units H — c — 1 are used. We present our results in Table 
4. Hence, one can see that there are degenerate eigenvalues of the spin partners within 
the attractive scalar and repulsive vector Yukawa potentials. For example, the Dirac spin 
doublet eigenstate partners are: {ls 1/2 , 0p 3/2 , 2p 3/2 , ld 5/2 , 3rf 5/2 , 2/ 7 / 2 ,4/ 7 / 2 , 3g 9/2 , 4/in /2 ), 
(2si /2 , 3p 3/2 , 0d 5/2 , 4d 5/2 , l/ 7/2 , 2g 9/2 , 3h 11/2 ), and (3si /2 , 4p 3/2 , 0/ 7/2 , lg 9/2 , 2h 11/2 , lp 1/2 , 0d 3/2 ). 
etc. 

Using Eq. (46), we can also calculate some relativistic particle E+ K and antiparticle 
E~ K pseudospin symmetric bound state energies with values of parameters a = 0.1 /to -1 , 
A — 1.0, M — 5.0 fm^ 1 and C ps = —5.0 /m _1 for various states with quantum numbers n 
and / (k < 0, k > 0) in the units h — c—1. We present our results in Table 5. One can see 
that there are degenerate eigenvalues of the pseudospin partners within the attractive scalar 
and repulsive vector Yukawa potentials. As an example, the Dirac pseudospin doublet 
eigenstate partners are: (lp 3/2 , 2si /2 , 2rf 5/2 ), (ld 5 / 2 , 2/ 7/2 , 0d 3/2 ), (lf 7/2 ,2g 9/2 ,0f 5/2 ,ld 3/2 ), 
(109/2, 2/in/2, 0^7/2, 1/5/2), and (1/111/2,0/19/2,1^7/2),- etc. When n = 1, there is a finite 
number of bound states where 1 < / < 25 (k — —25, k = 24) and when n = 2, 1 < / < 26 
(« = -26, k = 23) . 
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V. CONCLUSION AND OUTLOOK 

We have investigated the approximate bound state solutions of the Dirac equation for the 
screened Coulomb potential model with an arbitrary spin-orbit k- state under the conditions 
of the spin (pseudospin) symmetry limitation in the framework of the shortcut of the NU 
method including the usual approximation scheme to approximate the centrifugal (pseudo- 
centrifugal) barrier term. By setting S(r) (A(r)) to the spherically symmetric screened 
Coulomb potential model, we have derived the Dirac bound state energy spectra and asso- 
ciated two-component spinor wave functions for arbitrary spin-orbit K state that provides 
an approximate solution to the spin- and pseudo-spin symmetry. We have also discussed 
in detail how to choose the appropriate physical values for the spin constant C s from Eq. 
(25). In the presence of spin symmetry, we find that the appropriate value for C s falls in the 
range —20 jvrT x < C s < 8.8 fm~ l for positive energy part while jmT x < C s < 9.6 /m _1 
for the negative energy part. Thus, the allowed values for the C s in considering the whole 
spectrum (particle and antiparticle) falls in 4.9 fvrC 1 < C s < 8.8 jvrT x after neglecting the 
negative part of energy according to the requirement E n>K ^ —M. However, in the presence 
of pseudosymmetry, we find C ps value falls in the range —8 /m _1 < C ps < 20 fm~ 1 in the 
antiparticle energy spectrum demanding that E n>K ^ M. The Yukawa interaction between 
an electron and heavy nucleus appears to be less attractive as the value of a increases. The 
resulting solutions of the wave functions are being expressed in terms of the Jacobi polyno- 
mials (or hypergeometric functions). We have shown that the present spin symmetry can be 
easily reduced to the non-relativistic solution when one inserts S(r) = +V(r) (i.e., A(r) = 
or C s = 0). The non-relativistic limits of our solution are obtained by imposing appropriate 
changes of parameters k,(k +1) — *•/(/ + 1) in the spin symmetry limits. Furthermore, when 
a — > 0, our results can be reduced to the well-known bound state solutions for the Coulomb 
potential model. We must point out that the numerical calculations for eigenenergies of the 
Dirac states involved in Eqs. (28) and (46) are sensitive to the choice of the parameters C s , 
Cp S , a, A and M. The spin (pseudospin) limit Dirac energy spectrum computed in Table 3 
(Table 4) for arbitrarily chosen set of parameters are in the form of valence (hole) states. 
In order to remove the extra degeneracies in energy levels, we have solved Dirac- Yukawa 
problem in the presence of spin and pseudospin symmetry by adding a centrifugal-like term, 
i.e., V(r) = -Ae~ ar /r + D/r 2 . 



22 



Finally, the solution of the Dirac-Coulomb problem can be readily obtained from our 
solutions by setting a = 0. Hence, we can obtain expressions for the exact energy eigenvalues 
for the exact spin and pseudospin limitations. These exact solutions are identical to the ones 
found recently in Refs. [48,49]. 
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Appendix A: Resume of Parametric Generalization of NU Method 

We present the Nikiforov-Uvarov essential polynomials, root, eigenvalues and wave func- 
tions expressed in terms of the constants q (i — 1, 2, • • • , 13) and £^ (j= 0, 1, 2): 

(i) Constants: 

1 1 

c 4 = - (1 - ci) , c 5 = - (c 2 - 2c 3 ) , c 6 = cj + £ 2 , 

c 7 = 2c 4 c 5 - £i, c 8 = c 2 A + f , c 9 = c 3 (c 7 + c 3 c 8 ) + c 6 , 

,— 2 ^ 
cio = ci + 2c 4 + 2^ - 1 > -1, en = 1 - d - 2c 4 H Vc 9 > -1, 

c 3 
Cl2 = c 4 + v^ > 0, C13 = -c 4 H (v/Cg" - c 5 ) > 0, c 3 7^ 0. (Al) 

C3 

(ii) Polynomials: 

TT^) = C4+y/cs + C 5 X- (y/c^ + C 3 ^)x, (A2) 

fc = - (c 7 + 2c 3 c 8 ) - 2 V ^ 8 ^, (A3) 

r(x) = 1 + 2^ - [c 2 + 2 (^ + 03^ - c 5 )] x, (A4) 

r'(x) = -2c 3 - 2 (V^ + c 3v ^) < 0. (A5) 
(iii) Energy equation: 

c 2 n - (2n + 1) c 5 + (2n + 1) (^ + c 3y /c^) + n(n - l)c 3 + c 7 + 2c 3 c 8 + 2 V ^ 8 ^ = 0. (A6) 

(iv) Wave functions: 

p{x) =x Cl0 (l-c 3 x) Cn , (A7) 
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(j,{x) = X C12 (1 - C 3 X) C1S , C12 > 0, C13 > 0, (A8) 

y n (x) = Pi Cl0 ' Cll) (l - 2c 3 x), cio > -1, en > -1, (A9) 

Fn«(x) = A/; K a; Cl2 (l - c 3 x) Cl3 Pi cl0 ' c "Hl " 2c 3 x), 
= AT„ K x C12 (1 - c 3 x) Cl3 2 Fi (-n, 1 + ci + di + n; ci + 1; c 3 x) , x 6 [0, l/c 3 ] , c 3 ^ 0. (A10) 
where A/"„ K is a normalization constants. Further, P„ (x), ji> — 1,1/ > — 1 and x G [—1,1] 



are the Jacobi polynomials with 



P(^)(1-2 S ) 



fa + 1) 



n! 



jPi (— n, 1 + a + /3 + n; a + 1; s) 



(All) 



Appendix B: Calculations of the Normalization Constants 



The normalization constant, J\f n i can be determined in closed form. We start by using 
the relation between the hypergeometric function and the Jacobi polynomials (see formula 
(8.962.1) in [45]): 

1 — x 



Fi [ —n, n + v + /i + 1; 1/ + 1 



nl 



(* + l) 

to rewrite the wave functions in (32) as 



2 ) (iz + 1), 
r(n + i/ + l) 



-P^\x), 



r(i/ + i) ' 



(Bl) 



F nK (r) = ^n K ^2^+lf ~ 2Uiar{1 ~ e- 2a r +1 Pr? Ul ' 2K+1) (l ~ 2e- 2ar ). (B2) 

From the normalization condition f [u nj i(r)] dr = 1 and under the coordinate change 
x = 1 — 2e~ 2ar , the normalization constant in (B2) is given by 



AC 



-2 



1 



n!r(2i/i + 1) 
r(n + 2i/i + l) 



l2 .1 



-1 



2;/ 1 



2k+1 



1 — X \ / 1 + X \ I 1 + X 



[Pi 2 ^ 2K+1 \x)] 2 dx. 
(B3) 



The calculation of this integral can be done by writing 



l + x 



1 — X 
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and using the following two integrals (see formula (7.391.5) in [45]): 



(1 - xf- 1 (1 + xf [P^\x)] 2 dx = 2 v+ < i 



F(n + v + l)r(n + n + 1) 



n\vT(n + v + /j, + 1) 
which is valid for Re(z/) > and Re(/x) > —1 and (see formula (7.391.1) in [45]): 



(B4) 



(1 - xf (1 + xf [P { n v ^\x)] dx = 2 v+ » +l - 



(B5) 



_-, n\T(n + v + jj + l)(2n + z/ + /i + 1) 

which is valid for Re(u) > — 1, Re(/x) > — 1. Finally, we have carried out relativistic and 
non-relativistic normalization constants as 



K 

and 



r(2i/i + 1) 



i 



avAn + vi + k + 1) T(n + 1v x + l)F(n + 1v x + Ik + 2 



n + /t + l 



n!r(n + 2/t + 2) 



(B6) 



r(2 £n/ + 1) 

respectively. 



ae nl {n + I + e nl + 1) T(n + 2e nl + l)r(n + 2e n[ + 2/ + 2) 



(n + Z + 1} 



n!r(n + 2/ + 2) 



1/2 



(B7) 
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FIG. 1: Spin symmetry Dirac ground valence state energy level for (a) a particle and (b) anti- 
particle of mass M = 5.0 fm^ 1 as a function of the coupling constant A for several values of the 
range of the screened Coulomb potential. 



FIG. 2: Dirac ground valence state energy level for (a) a particle and (b) anti-particle of mass 
M = 5.0 /to -1 as a function of the spin symmetry constant C s for several values of the spin-orbit 



K. 



FIG. 4: Pseudospin symmetry Dirac ground hole state energy level for (a) a particle and (b) anti- 
particle of mass M = 5.0 fm^ 1 as a function of the coupling constant A for several values of the 
range of the screened Coulomb potential. 



FIG. 5: Dirac ground hole state energy level for (a) a particle and (b) anti-particle of mass M = 5.0 
fm^ 1 as a function of the pseudospin symmetry constant C ps for several values of the spin-orbit 

K. 



FIG. 3: Spin symmetry Dirac upper and lower spinor wave functions of ground Opi/2 and first 
excited lpi/2 states of (a) particle and (b) antiparticle with n = 1. 
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FIG. 6: Pseudospin symmetry Dirac upper and lower spinor wave functions of ground 0d 3 / 2 and 
first excited W3/2 states of (a) particle and (b) antiparticle with k = 2. 



TABLE I: Specific values of the constants in Dirac- Yukawa problem. 



Spin symmetry 


case 




Pseudospin symmetry case 


Cl = l 






Cl = 


1 


C2 = l 






C2 = 


1 


C3 = l 






C3 = 


1 


c 4 = 






C 4 = 





C 5 = -^ 






C5 = 


1 
2 


c 6 = i [1 + 4 (z/ : 


i+w 


■l)] 


C6 = 


i[l + 4(^+ W2 )] 


C7 = — 2v\ — oj\ 


+ k( 


K + l) 


Cl = 


— 2v\ — 0J2 + k {k — 1) 


c 8 = v\ 






C8 = 


"2 


c 9 = i (2k + l) 2 






C9 = 


I (2k - l) 2 


cio = 2i/i 






ClO = 


= 2z/ 2 


Cn = 2k + 1 






Cll = 


= 2k- 1 


C12 = ^1 






C12 = 


= ^2 


C13 = K + 1 






Cl3 = 


= K 


£ 2 = v\ + Wl 






e 2 = 


: v\ + w 2 


^ = 2i/? + wi- 


K (K 


+ 1) 


£1 = 


: 2z^2 + ^2 — K (k — 1) 


*o = "? 






e = 


v\ 



30 



TABLE II: Spin symmetric Dirac- Yukawa ground valence energy spectrum (in units of fm 1 ) of 
the particle and antiparticle as a function of C s for various values of k. 

C s Eq,k=1 ^0,k=3 ^0,k=5 ^0,k=7 ^0,k=9 ^s ^o,k=1 ^0,k=3 ^0,k=5 

-20 3.328 4.632 4.880 4.962 4.992 -20 -24.999 -24.995 -24.988 

-15 3.622 4.707 4.913 4.977 4.998 -15 -19.998 -19.994 -19.986 

-10 3.916 4.783 4.943 4.990 5.000 -10 -14.998 -14.992 -14.982 

-5 4.209 4.857 4.972 4.999 4.996 -5 -9.997 -9.989 -9.976 

4.502 4.929 4.995 4.997 4.975 

5 4.792 4.990 4.993 4.951 4.883 1 

6 4.849 4.998 4.982 4.921 4.829 4 

7 4.905 5.000 4.958 4.865 4.726 4.8 

8 4.957 4.988 4.897 4.722 4.363 4.9 

8.5 4.980 4.968 4.818 - - 5 0.007856 0.032 0.072 

8.8 4.992 4.942 4.689 6 1.010 1.039 1.090 



4.996 


-4.984 


-4.964 


3.996 


-3.982 


-3.960 


0.993 


-0.974 


-0.940 


0.192 


-0.170 


-0.131 


0.092 


-0.069 


-0.030 



9 4.998 


4.910 


— 


— 


— 


7 


2.013 


2.053 


2.121 


9.5 4.987 


— 


— 


— 


— 


8 


3.019 


3.079 


3.188 


10 - 


— 


— 


— 


— 


9 


4.038 


4.173 


— 


10.5 - 


— 


— 


— 


— 


9.5 


4.577 


— 


— 


10.6 5.498 


— 


- 


— 


— 


9.6 


4.702 


— 


— 


10.7 5.623 


— 


- 


— 


— 


9.7 


— 


— 


— 


10.8 5.737 


— 


— 


- 


— 


10.5 


— 


— 


— 


10.9 5.846 


— 


- 


— 


- 


10.6 5.232 


— 


— 


11 5.953 


5.754 


— 


— 


— 


11 


5.200 


5.360 


— 


11.1 6.058 


5.899 


— 


— 


— 


11.3 5.205 


5.272 


— 


11.3 6.266 


6.147 


— 


— 


— 


11.4 5.208 


5.258 


5.511 


11.4 6.368 


6.262 


5.998 


— 


— 


11.5 5.212 


5.247 


5.450 


12 6.979 


6.912 


6.786 


6.562 


— 


12 


5.233 


5.218 


5.327 


13 7.986 


7.944 


7.870 


7.759 


7.600 


15 


5.396 


5.208 


5.208 


14 8.990 


8.959 


8.905 


8.828 


8.724 


20 


5.686 


5.268 


5.205 


15 9.992 


9.967 


9.925 


9.866 


9.787 










20 14.996 


14.984 


14.964 


14.935 


14.898 
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TABLE III: Nonrelativistic bound state energies of the Yukawa potential in units of % = m = 1. 
For comparison, we set A = \/2 and a = gA. 



State nig NU (Present) [13] (Numerical) [23] (AIM) [21] (SUSY) [16] 



[17] 



Is 


0.002 


-0.996004 


-0.9960 


-0.996006 


-0.99601 


-0.99601 


-0.9960 




0.005 


-0.990025 


-0.9900 


-0.990037 


-0.99004 


-0.99004 


- 




0.01 


-0.9801 


-0.9801 


-0.980149 


-0.98015 


-0.98015 


-0.9801 




0.02 


-0.9604 


-0.9606 


-0.960592 


-0.96059 


-0.96059 


-0.9606 




0.025 


-0.950625 


-0.9509 


-0.950922 


-0.95092 


-0.95092 


- 




0.05 


-0.9025 


-0.9036 


-0.903632 


-0.90363 


-0.90363 


-0.9036 


2s 


1 0.002 


-0.246016 


-0.2460 


-0.246023 


-0.24602 


-0.24602 


-0.2460 




0.005 


-0.2401 


-0.2401 


-0.240148 


-0.24015 


-0.24015 


- 




0.01 


-0.2304 


-0.2306 


-0.230586 


-0.23059 


-0.23058 


-0.2306 




0.02 


-0.2116 


-0.21230 


-0.212296 


-0.21230 


-0.21229 


-0.2124 




0.025 


-0.2025 


-0.2036 


-0.203551 


-0.20355 


-0.20352 


- 




0.05 


-0.16 


-0.1635 


-0.163542 


-0.16351 


-0.16325 


-0.1650 


2p 


1 0.002 


-0.246016 


-0.2460 


-0.246019 


-0.24602 


-0.24602 


-0.2460 




0.005 


-0.2401 


-0.2401 


-0.240123 


-0.24012 


-0.2412 


- 




0.01 


-0.2304 


-0.2305 


-0.230490 


-0.23049 


-0.23049 


-0.2305 




0.02 


-0.2116 


-0.2119 


-0.211926 


-0.21192 


-0.21193 


-0.2120 




0.025 


-0.2025 


-0.2030 


-0.202984 


-0.20299 


-0.20299 


- 




0.05 


-0.16 


-0.1615 


-0.161480 


-0.16144 


-0.16155 


-0.1625 


3p 


1 1 0.002 


-0.107147 


-0.1072 


-0.107160 


-0.10716 


-0.10716 


-0.1072 




0.005 


-0.101336 


-0.1014 


-0.101416 


-0.10142 


-0.10142 


- 




0.01 


-0.092011 


-0.09231 


-0.092306 


-0.09231 


-0.09236 


-0.09236 




0.02 


-0.074711 


-0.07570 


-0.075704 


-0.07570 


-0.07563 


-0.07611 




0.025 


-0.066736 


-0.06816 


-0.068157 


-0.06814 


-0.06799 


- 




0.05 


-0.033611 


-0.03712 


-0.037115 


-0.03739 


-0.03486 


-0.04236 


3d 


2 0.002 


-0.107147 


-0.1072 


-0.107152 


-0.10715 


-0.10715 


-0.1072 




0.005 


-0.101336 


-0.1014 


-0.101368 


-0.1014 


-0.10137 


- 




0.01 


-0.092011 


-0.09212 


-0.092122 


-0.09212 


-0.09212 


-0.09216 




0.02 


-0.074711 


-0.07503 


-0.075030 


-0.07502 


-0.07504 


-0.07531 




0.025 


-0.066736 


-0.06715 < 


jo -0.067146 


-0.06713 


-0.06718 






TABLE IV: Dirac valence bound state energies (in units of fm 1 ) of the spin symmetry Yukawa 
problem for various values of n and I. 



1 n 


, K < 


(I, j) 


E n ,K<o E nK<0 n 


, K > 


{I, j) 


^n,K>0 ^n,K>0 





-1 


0*1/2 


- - 






1 


-1 


ls l/2 


-0.098076 4.05808 






2 


-1 


2si/ 2 


-0.0922956 4.78641 






3 


-1 


3si/ 2 


-0.0826327 4.94209 






4 


-1 


4si/ 2 


-0.0690436 4.98904 






1 


-2 


0^3/2 


-0.098076 4.05808 


1 


QPl/2 


-0.0922956 4.78641 


1 1 


-2 


1P3/2 


1 


1 


lPl/2 


-0.0826327 4.94209 


1 2 


-2 


2p 3 /2 


-0.098076 4.05808 2 


1 


2^1/2 


-0.0690436 4.98904 


1 3 


-2 


3P3/2 


-0.0922956 4.78641 3 


1 


3P1/2 


-0.0514655 4.99998 


1 4 


-2 


4p 3 /2 


-0.0826327 4.94209 4 


1 


4pi/ 2 


-0.0298154 4.99395 


2 


-3 


04/2 


-0.0922956 4.78641 


2 


04/2 


-0.0826327 4.94209 


2 1 


-3 


14/2 


-0.098076 4.05808 1 


2 


14/2 


-0.0690436 4.98904 


2 2 


-3 


2 4/2 


2 


2 


2 4/2 


-0.0514655 4.99998 


2 3 


-3 


34/2 


-0.098076 4.05808 3 


2 


34/2 


-0.0298154 4.99395 


2 4 


-3 


44/2 


-0.0922956 4.78641 4 


2 


44/2 


-0.0039874 4.97759 


3 


-4 


0/7/2 


-0.0826327 4.94209 


3 


0/5/2 


-0.0690436 4.98904 


3 1 


-4 


1/7/2 


-0.0922956 4.78641 1 


3 


1/5/2 


-0.0514655 4.99998 


3 2 


-4 


2/7/2 


-0.098076 4.05808 2 


3 


2/5/2 


-0.0298154 4.99395 


3 3 


-4 


3/7/2 


3 


3 


3/5/2 


-0.0039874 4.97759 


3 4 


-4 


4/7/2 


-0.098076 4.05808 4 


3 


4/5/2 


-0.0261492 4.95362 


4 


-5 


0^9/2 


-0.0690436 4.98904 


4 


0#7/2 


-0.0514655 4.99998 


4 1 


-5 


109/2 


-0.0826327 4.94209 1 


4 


107/2 


-0.0298154 4.99395 


4 2 


-5 


2 09/2 


-0.0922956 4.78641 2 


4 


2 07/2 


-0.0039874 4.97759 


4 3 


-5 


3ff9/2 


-0.098076 4.05808 3 


4 


307/2 


-0.0261492 4.95362 


4 4 


-5 


409/2 


4 


4 


4^7/2 


-0.0607542 4.92325 


5 


-6 


0/tll/2 


-0.0514655 4.99998 


5 


O/lg/2 


-0.0298154 4.99395 


5 1 


-6 


1^11/2 


-0.0690436 4.98904 1 


5 


14/2 


-0.0039874 4.97759 


5 2 


-6 


2/i n /2 


-0.0826327 4.94209 2 


5 


2/l 9 /2 


-0.0261492 4.95362 


5 3 


-6 


3/111/2 


-0.0922956 4.78§^1 3 


5 


34/2 


-0.0607542 4.92325 



TABLE V: Dirac hole bound state energies (in units of fm' 
problem for various values of n and I. 



x ) of the pseudospin symmetry Yukawa 



I n 


, K < 


{I, J) 


^n,K<0 


Kn<o ri-1, 


k >0 (1 + 2, 


i + 1) E n, K >0 


- C/ n,K>0 


1 1 


-1 


!si/2 


- 


0,2 


0^/2 


-4.603587 


-0.00817777 


2 1 


-2 


lj»3/2 


-3.917958 


-0.00204188 0,3 


0/5/2 


-4.749130 


-0.01843870 


3 1 


-3 


14/2 


-4.603587 


-0.00817777 0,4 


Os-7/2 


-4.791740 


-0.03287710 


4 1 


-4 


1/7/2 


-4.749130 


-0.01843870 0,5 


O/19/2 


-4.799920 


-0.05156860 


5 1 


-5 


159/2 


-4.791740 


-0.03287710 0,6 




-4.791590 


-0.07461340 


6 1 


-6 


1^11/2 


-4.799920 


-0.05156860 0,7 




-4.772990 


-0.102140 


7 1 


-7 




-4.791590 


-0.07461340 0,8 




-4.746630 


-0.134308 


8 1 


-8 




-4.772990 


-0.102140 0,9 




-4.713610 


-0.171314 


9 1 


-9 




-4.746630 


-0.134308 0,10 




-4.674380 


-0.213399 


10 1 


-10 




-4.713610 


-0.171314 0,11 




-4.629040 


-0.260854 


15 1 


-15 




-4.454290 


-0.439464 0, 16 




-4.300520 


-0.594696 


20 1 


-20 




-3.993200 


-0.903412 0,21 




-3.708850 


-1.188370 


25 1 


-25 




-2.938620 


-1.959250 0,24 




-2.938620 


-1.959250 


1 2 


-1 


2«l/2 


-3.917958 


-0.00204188 1,2 


14/2 


-4.749130 


-0.01843870 


2 2 


-2 


2^3/2 


— 


1,3 


1/5/2 


-4.791740 


-0.03287710 


3 2 


-3 


24/2 


-3.917958 


-0.00204188 1,4 


107/2 


-4.799920 


-0.05156860 


4 2 


-4 


2/7/2 


-4.603587 


-0.00817777 1, 5 


14/2 


-4.791590 


-0.07461340 


5 2 


-5 


2^9/2 


-4.749130 


-0.01843870 1,6 




-4.772990 


-0.102140 


6 2 


-6 


2/&11/2 


-4.791740 


-0.03287710 1,7 




-4.746630 


-0.134308 


7 2 


-7 




-4.799920 


-0.05156860 1,8 




-4.713610 


-0.171314 


8 2 


-8 




-4.791590 


-0.07461340 1,9 




-4.674380 


-0.213399 


9 2 


-9 




-4.772990 


-0.102140 1,10 




-4.629040 


-0.260854 


10 2 


-10 




-4.746630 


-0.134308 1,11 




-4.577470 


-0.314039 


15 2 


-15 




-4.519370 


-0.373394 2, 16 




-4.209880 


-0.685887 


20 2 


-20 




-4.108150 


-0.788067 2,21 




-3.525210 


-1.372260 


25 2 


-25 




-3.291060 


-1.606630 2,23 




-2.938620 


-1.95925 


26 2 


-26 




-2.938620 


-1.959250 
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